Controlled-NOT logic with nonresonant Josephson phase qubits 
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We establish theoretical bounds on qubit detuning for high fidelity controlled-NOT logic gate 
implementations with weakly coupled Josephson phase qubits. It is found that the value of qubit 
detuning during the entangling pulses must not exceed 2g for two-step, and g for single-step control 
sequences, where g is the relevant coupling constant. 
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I. INTRODUCTION 

In ourprevious work on steering with Josephson phase 
qubits [l| we found two peculiar single-step controlled- 
NOT (CNOT) implementations involving off-resonance 
qubits detuned by the amount smaller than the char- 
acteristic coupling constant. Does the detuning always 
have to be so small? What happens if it gets larger? 

These questions seem to be important on several 
counts. 

Firstly the majority of entangling gate designs pro- 
posed in the literature assume resonant qubits @, 0, 0, 
S, S, 0, B 1- (However, see Refs. [13, El for notable 
exceptions.) This is hardly surprising since in the rotat- 
ing wave approximation (RWA) typically used to analyze 
superconducting qubits the resonant condition leads to 
relatively simple and easily solvable Hamiltonians con- 
taining no local crj? terms. This works well when the sys- 
tem consists of only two qubits. However, it is reason- 
able to expect that when thousands of such qubits are 
assembled into an integrated circuit, maintaining them 
on resonance will become a difficult task. How then can 
we be sure that an architecture involving detuned qubits 
is able to reliably generate universal gates, such as, for 
example, a CNOT gate? 

Secondly, some qubits may be fabricated with various 
defects preventing them from being tuned to resonance 
exactly. Would they still be usable? And what if for some 
applications it becomes advantageous (or even necessary) 
to use qubits with intrinsically different level splittings? 

Thirdly, we should also keep in mind that detuning is 
routinely used in actual experiments to "decouple" the 
qubits in order to perform local (that is, nonentangling) 
operations @, [H, [lj, [3| . It would be useful if, after such 
decoupling is performed, we do not have to worry about 
bringing the qubits back to exact resonance when doing 
subsequent entangling operations. 

Additionally, from the purely theoretical viewpoint, 
decoupling provides a useful limit against which to check 
our calculations. If at larger detuning the interaction is 
expected to lose its entangling properties, we must be 



able to predict a crossover into the regime when the gate 
fidelity (of the CNOT, in our case) starts to deteriorate. 
Thus, the primary goal of this paper will be to establish 
the exact conditions under which such crossover occurs 
for previously proposed CNOT implementations involv- 
ing superconducting qubits. 



The paper is organized as follows: 



In Section 111X1 we introduce the Hamiltonian for capac- 
itively and inductively coupled Josephson phase qubits. 
Since there is an infinite number of possible local drives 
and rotating frames to choose from, it will be necessary 
to limit our discussion to situations that are simple and 
yet powerful enough to provide some new and interest- 
ing physical insights. Our specific choices for Rabi pulses 
and rotating frames will be made in Eqs. ((3]), ((4]), and 
of that Section. 



In Section IIVI we generalize to finite detuning the fa- 
miliar two-step CNOT sequence involving a local 7r-pulse 
sandwiched between two entanling operations. The ex- 
act bound on detuning that guarantees generation of the 
perfect (in the RWA) controlled-NOT logic gate will then 
be given in Eq. ([?!]) . 



In Section |V] we numerically solve the single-step case. 
We will see that in that case the restriction on detuning 
for perfect CNOT generation is somewhat stronger than 
in the two-step case (Eq. (f!?T|) ). 



In Section I VII we present optimal results for approx- 
imate CNOT gates at detuning larger than maximally 
allowed. We will see how the Makhlin invariants and 
fidelity of the optimized gates deviate from their ideal 
CNOT values. We will also simulate the Weyl cham- 
ber steering trajectories corresponding to such optimized 
gates. 



'Electronic address: ag@physast.uga.edu 



We conclude in Section I VIII with a brief summary of 
our results. 
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II. NOTATION 

In what follows, we will use the notation that is con- 
venient for Lie algebraic manipulations, 



Xi. 



-a\, Z k = -o%, (fc = l,2), 



XX = -ofo* YY = -a%af, ZZ = -a z 2 af, 
XY= l -***y, YX= i -ataf. 



(1) 



Notice, that [XX, YY] = [XY,YX] = 0, and ZZ com- 
mutes with each of the XX, YY, XY, YX operators. 

Why is this notation convenient? 

Consider the following transformation (called "going 
to a rotated frame") on the Lie algebra su(4) of the two- 
qubit system: 



XX 



- 0Zl XXe eZl = XXcos9 + XYsmi 



(2) 



This transformation can be nicely interpreted as a rota- 
tion of vector XX by an angle 9 in the real vector space 
spanned by the generators of the group SU(4). This is 
how the continuous group acts on its Lie algebra. Math- 
ematicians call it the adjoint representation. The alge- 
bra plays the role of the representation space for its own 
group. 

Since we do work at the level of algebra, and not at 
the level of the Hilbert space when discussing equivalence 
classes of gates, this is a very convenient notation. It sim- 
plifies things. Also notice how naturally the periodicity 
of 2ir appears at this level of description. 

Now, if we were to write the same transformation in 
terms of the Pauli matrices, we would have to remember 
to put in the imaginary unit i and the factors of 1/2 in 
the exponents on the left hand side of Eq. ([2]), 



iOa 1 12 ( _x x\ 



i9<rl/2 _ 



(crf(7i)cos0 + {alal) sin 9, 



breaking its beautiful symmetry. 



III. THE HAMILTONIAN 

When restricted to the computational subspace, the 
Hamiltonian for two coupled Josephson phase qubits, one 
of which is driven by a resonant rf pulse, is given by 

iH(t) = -ojZi - (u + 5)Z 2 

+2fii cos (tot) Xj. + 2(gYY + gZZ), (3) 

where g,g<^uj are the coupling constants, to is the level 
splitting of the first qubit, \5\ <C uj is the detuning, and 
Oi is the corresponding Rabi frequency. Here we assume 
that g ~ Qi, which differs from the condition g -C il\ <C 
to adopted in Ref. [l(|. For realistic systems, u> « 10 
GHz, g « 10 MHz. For capacitive coupling, g = 0; for 
inductive coupling, g < O.lg. 



In the doubly rotating frame defined by 
e lHot (...)e- tHot , with iH = -u{Z x + Z 2 ), after 
averaging over fast oscillations, the system Hamiltonian 
is time- independent, 



RWA 



= -SZ 2 + QiXi + iU + gZZ, 



with 



iU = g(XX + YY) = ig 



/0 0\ 

10 

10 

\0 0/ 



(4) 



(5) 



Alternatively, to perform a useful consistency 
check, we consider another rotating frame defined by 

e i6ot (. . . )e~ lAot , with iH = -uZ x - (w + 6)Z 2 . In this 
frame, the RWA Hamiltonian is 



iH BWA (t) = £l l X 1 +iH(t)+gZZ, 



(6) 



where now we have a slowly varying interaction term 
given by 



M{i) 



g [(XX + YY) cos(5t) + (YX - XY) sin(St)] 

/0 0^ 

e' lSt 
e +lSt 



(7) 



\0 0/ 



The central block of this matrix has the form of a rotat- 
ing drive for a spin-1/2 system for which the analytical 
solution is well known [191 ]. This observation will prove 
helpful for calculations in Section [TV] 

We will now show how these two RWA Hamiltonians 
lead to locally equivalent CNOT implementations. 



IV. TWO-STEP CNOT 

The well-known two-step CNOT implementation for 
resonant (S — 0) qubits has control sequence [HI, [l6[ 

CNOT (2) = e l(7r/4) i?p OSt [U(t (2) )e-^U(t (2) )] R pic , (8) 

with ip) — 7r /4ff, where 

U(t (2) ) = e -H2)(9(XX+YY)+gZZ) 



( 



0\ 



l/\/2 -i/y/2 
-i/y/2 1/V2 
\0 1/ 



, (9) 



and i? P ost,pre are some local rotations. We will choose 



R 



post 



pre 



_ -(7r/2)Z 2 +(ir/2)(X 2 +X 1 ) 



Of particular importance to us is the entangling part 
U(t)e~ 7rXl U(t) that determines the local equivalence 
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class of the full gate. In our case the local class is 
controllcd-NOT whose canonical representative in the 
computational basis is defined to be 



CNOT ee 



'1 0\ 

,!, 1 € U ( 4 )' det(CNOT) = -1. 
vO 1 0/ 



(11) 

We ask the following question: If \S\ > (detuned 
qubits), can we still use the sequence in Eq. (jSJ) to gener- 
ate a CNOT gate, possibly with different gate time and 
different pre- and post-rotations? The answer to this 
question turns out to be "Yes" , provided S is restricted 
in a certain way. 

Setting fix = in Eqs. (|3]) and © gives the time 
evolution operator 



U(t) = e 



in frame 1, and 



fe m l 2 \ 

u —iv 

— iv u* 

\ e ~ lSt l 2 ) 



(12) 



U{t) = e 



-tgZZ 



A o oo^ 

ue~ lSt / 2 -ive- i5t ' 2 

-ive lSt / 2 u*e iSt / 2 

\0 1/ 



(13) 



in frame 2, where 



/ JS 2 + Ag 2 \ iS 
u = cos — 1 -I , sm 

V 2 j v^ + V 



Vs~ 2 



(14) 



2.9 



^S 2 + V 



: sm — i 



(15) 



In both frames the Makhlin invariants [17[ of 
U{t)e-* x ^U{t) are 



V 



S 2 



/ 



(16) 



and 



G 2 = (35 4 + 85 2 .g 2 l + 2cos(tv/j 2 +4 5 2 ) 

+16 5 4 [2 + cos(2iv/^TV)] )/ (S 2 + Ag 2 )' , 

(17) 

and are independent of the ZZ coupling. This shows that 
for any t the resulting gates are represented by the same 
point on the (XX, Y"Y")-plane of the Weyl chamber (see 
[lL Hq for discussion). Since CNOT class corresponds to 
G\ = 0, G2 = 1, we get 



t(£\ — 



7r — arccos 



(S 2 /4g 2 



(18) 



^/S 2 + 4<? 2 

with the limit t( 2 ) ^ ^/Ag trivially recovered for vanish- 
ing detuning. 

For example, for g = 0, 5 = l.OOg, we get i( 2 ) = 
1.03837r/4(7. The corresponding CNOT gate is given by 
Eq. ([51), where now 



U(t {2 



/0.9180 + 0.3965i \ 

0.6124 + 0.3536i -0.7071i 

-0.7071* 0.6124 - 0.3536i 

V 0.9180- 0.3965i/ 



R 



post 



e -(7r/2)y 2e -(7r/2)(a 2 Z2+ai^i) 



Rpre 



e -(ir/2)((l+a 2 )Z 2 +a 1 Z 1 ) e +(Tr/2)(X 2 +X 1 ) 



(19) 



with a2 = 0.5929, ot\ = 0.2596 in frame 1, and 



U(t { 



2)j 



a o o\ 

0.7024 + 0.0817i -0.2804- 0.6491? 

0.2804- 0.6491i 0.7024 - 0.0817i 

v0 1/ 



the perfect in the rotating wave approximation canon- 
ical controlled-NOT gate have been found numerically. 
Notice that the two-step control sequence works only for 



R 



post 



e -(7r/2)Y 2e -W2)/3Z 2 



d _ -(■n/2)(l+fi)Z 2 +(^/2){X 2 +X 1 ) 
1 *-pre — c c j 



(20) 



with ft = -0.1858, = 0.3333 in frame 2. Here, the 
pre- and post-rotations (cf. Eq. (flTj|) ) needed to generate 



|<5| < 2 ff < u). 



SINGLE-STEP CNOT 



(21) 



Here, for simplicity, we limit our discussion to the RWA 
Hamiltonian given in Eq. (|U) with g — (capacitive 
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coupling). [In this case, a single Rabi term suffices to 
implement a CNOT gate. When j ^ 0, an additional 
local Rabi drive ~ Q2X2 must be applied to the second 
qubit.] The CNOT sequence is then 



CNOT (1) = e^^R post U(t {1) )R v 



(22) 



where = tt/2 5 , = g^(4n) 2 - 1, n = 1, 2, 3, . . . , 



/ 1 -i\ 

(-1)" 1 -? 

y/2 -i 1 

\-i 1 



(23) 



and 



R , - p -(T/2)y2 D _ -(n/2)Z 2 +(tt/2)(X 2 -X 1 ) 

. (24) 

In order to generalize this single-step implementation to 
finite detuning, we optimize the gate parameters tn\ and 
fix of 



= p -t (1) [-SZ 2 +Q 1 X 1 +g(XX+YY)] 



(25) 



using Nelder-Mead simplex direct search with bound con- 
straints for the minimum of the distance from the CNOT 
class defined by 

d 2 (fti,i(i)) := |Gi(ni,t ( i))-Gi(C3NOT)| 2 

+ |G 2 (fii, i (1) )-G 2 (CNOT)| 2 , (26) 



with Gi ) 2(fii,t(i)) being the Makhlin invariants of 
{/(£(!)). It is important to keep in mind that the dis- 
tance function introduced in Eq. ()26|1 is not a measure 
of the gate fidelity. Infinitely many gates — all differing 
from each other by arbitrary local rotations — may have 
the same value of d 2 (fli, i(i))- Once the entangling part 
f/(t(i)) is found to have d 2 (fli, tm) — 0, it can then be 
made into the canonical CNOT gate by additional local 
rotations. 

Actual experiments motivate this choice of the distance 
function. It is generally believed that doing local rota- 
tions is easy, but performing entanglement is difficult. 
Thus, if by using experimentally available interaction and 
the local controls we can somehow steer the system into 
the "right" equivalence class, then making the actual tar- 
get gate will be relatively straightforward. 

Performing the optimization we find that in order to 
generate the exact CNOT in the single-step case the de- 
tuning must be restricted by 



\S\ <g<Luj. 



(27) 



For example, for maximally allowed 6 — g, we get 
t (1 ) = 1.27537r/20, fii = 3.7781g. The resulting gate 
is given by Eq. (f2"2")l with 



-0.2553- 0.4300?' 0.4821 - 0.1324?' -0.4821 + 0.1324?' 



0.5001?' 



U(t (1 



R 



post 



0.4821 -0.1324? -0.0001 + 0.5001?' 0.5001? 0.4821 + 0.1324?' 

-0.4821 + 0.1324?' 0.5001? -0.0001-0.5001? 0.4821 + 0.1324?' 

0.5001? 0.4821 + 0.1324i 0.4821 + 0.1324?' -0.2553 + 0.4300?/ 

{■n/2)(a 2 Z 2 +a 1 Z 1 ) R _ p - {tt /2)((l+a 2 )Z 2 +a 1 Z 1 ) p + /2)(X 2 ~(l+~/ 1 )X 1 ) 



pre 



(28) 



where a 2 = 0.8294, a x = -0.1705, 71 = -0.9998. 

To visualize how this gate is reached in the course 
of the unitary evolution we simulate the Weyl chamber 
steering trajectory c(t) = {c\ (t),c%(t), cs(t)) for the above 
mentioned values of the gate parameters. The goal here 
is to establish a correspondence 16], 

U(t) ~ e -ci(t)XX-C2(t)YY-c 3 (t)ZZ ^ , 2 9) 

between the physical gate U(t) and the unphysical ma- 
trix exponential that formally resides in the same local 
equivalence class as U (t) . The time-dependent vector c(t) 
then represents the dynamically generated class at every 
moment of system's evolution. Figures [T] and [2] show how 
the CNOT class that has c = (7r/2, 0,0) is generated in 
our single-step example with 5 = g. 



For other values of qubit detuning, the gate parameters 
required to implement perfect two-step and single-step 
CNOT logic are listed in Table H 



VI. WHAT HAPPENS AT LARGER 
DETUNING? 

When detuning exceeds the limits set in Eqs. (|2"Tj) and 
(|27[) . our implementations begin to deviate from their 
perfect controlled-NOT form. To visualize how that hap- 
pens, we again search for the gate parameters that pro- 
duce the gates that belong to the local equivalence classes 
closest to CNOT in the sense of the distance function de- 
fined in Eq. (|26[) . The corresponding results for single- 
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0.2 0.4 0.6 0.0 1 1.2 

t, units of n/2g 



FIG. 1: Time dependence of ci (top curve) and ci (bot- 
tom curve) for single-step CNOT implementation with capac- 
itively coupled Josephson phase qubits at maximal detuning, 
S = g. Here, cs — at all times. Dashed curves represent the 
resonant case, 8 = 0. 




FIG. 2: Single-step steering trajectory corresponding to Fig. 
[T] that generates CNOT class at maximal detuning, S = g. 
The steering parameters are given in units of tt/2. Dashed 
curve is for 5 = 0. 



TABLE I: Generation of ideal controlled-NOT logic (Gi = 0, 
G2 = 1) using capacitively coupled Josephson phase qubits at 
finite detuning. In the two-step case, \6\ < 2g <C u), fii = 0, 
and t(2) = (7r/4g)T(2), as given in Eq. (|18p . In the single-step 
case, \S\ < g <C w. The optimal values of fii w g\/l5 and 
t(x\ = (n /2g)T( 1 ) have been found numerically using Nelder- 
Mead simplex direct search with bound constraints. 



1*1/9 


^(2) 




T(i) 




0.00 


1.0000 




1.0000 


3.8730 


0.10 


1.0003 




1.0009 


3.8724 


0.20 


1.0014 




1.0037 


3.8707 


0.30 


1.0031 




1.0085 


3.8679 


0.40 


1.0056 




1.0155 


3.8638 


0.50 


1.0088 




1.0253 


3.8583 


0.60 


1.0128 




1.0386 


3.8513 


0.70 


1.0177 




1.0568 


3.8422 


0.80 


1.0235 




1.0827 


3.8303 


0.90 


1.0303 




1.1245 


3.8132 


1.00 


1.0383 




1.2753 


3.7781 


1.10 


1.0476 








1.20 


1.0585 








1.30 


1.0713 








1.40 


1.0863 








1.50 


1.1042 








1.60 


1.1261 








1.70 


1.1536 








1.80 


1.1901 








1.90 


1.2445 








2.00 


1.4142 








TABLE II: 


Optimized values of the gate parameters for the 


gates closest to the single-step CNOT in the sense of Eq. (|26p 


at \S\ > g. 


Here, tm = 




Notice how the Makhlin 


invariants deviate from their ideal CNOT values for larger S. 


\S\/g 


T M 


Oi/a 


Gi 


G 2 


1.00 


1.2753 


3.7781 


0.0000 


1.0000 


1.10 


1.2330 


3.7470 


0.0030 


0.9994 


1.20 


1.1945 


3.7323 


0.0106 


0.9978 


1.30 


1.1590 


3.7250 


0.0214 


0.9955 


1.40 


1.1262 


3.7203 


0.0340 


0.9927 


1.50 


1.0961 


3.7152 


0.0476 


0.9898 


1.60 


1.0686 


3.7074 


0.0614 


0.9867 


1.70 


1.0438 


3.6952 


0.0749 


0.9837 


1.80 


1.0216 


3.6772 


0.0879 


0.9808 


1.90 


1.0019 


3.6519 


0.1003 


0.9780 


2.00 


0.9849 


3.6179 


0.1118 


0.9754 



pulses (not shown) results in the optimized gate, 




0.9866 -0.1122? 0.0258i 0.1158 \ 

-0.1122i 0.9866 -0.1158 -0.0258i | 

-0.1186 0.0009i 0.1122i 0.9866 

-0.0009i 0.1186 0.9866 0.1122a / 



(30) 

whose intrinsic fidelity with respect to the canonical 
controlled-NOT gate is 

T = yjl-ti [(U opt - CNOT) t (U opt - CNOT)] 
= 0.9448. (31) 



step implementation are given in Table HT1 

Figures [3] and 2] show the simulated single-step steering 
trajectory generated by the RWA Hamiltonian given in 
Eqs. (HI) and (|2"5"|) at 6 = 1.5g. Direct search for local 



VII. CONCLUSION 

In summary, we have demonstrated that the CNOT 
pulse sequences previously proposed for resonant Joseph- 
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son phase qubits may still be used at finite detuning. To 
achieve high fidelity of the resulting gate the value of the 
detuning during the entangling operations should not be 
greater than 2g in the two-step implementation and g in 
the single-step implementation, respectively. 
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FIG. 3: Time dependence of ci (top curve) and ci (bot- 
tom curve) for single-step implementation generating the class 
closest to CNOT in terms of Eq. (TJSI at S = 1.5g. Here, 
C3 = at all times. Dashed curves represent the resonant 
case. Time is measured in units of Tt/2g. The closest class is 
reached at t(j^ = 1.0961. 
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FIG. 4: Single-step steering trajectory corresponding to Fig. 
[3] at detuning S — 1.5g. The steering parameters are given in 
units of 7r/2. Dashed curve is for 5 = 0. 
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